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[- — Abstract 

Using a duality approach we give an example of modeling nonlocal and mixed condensates 
r— i whose behavior mimics that of QCD. Quantitatively, the estimated value of the condensate 

parameter TOq is approximately 0.6 — 0.7 GeV 2 . 

1 Introduction 

The determination of the condensates of the QCD vacuum is a very important issue in phe- 
nomenology of strong interactions [lj. As known, it can only be done in a non-perturbative 
formulation of the theory. At present, not much is known about the condensates from first prin- 
ciples. Some attempts exist to determine them on the lattice or in the models of the instanton 

qq vacuum. In practice, however, the QCD sum rules still remain a basic tool for doing soQ Thus, 

there is a strong need for new alternative approaches to the problem. 

The AdS/CFT correspondence [3] has opened new avenues for studying strongly coupled 
• • gauge theories. Although the original proposal was for conformal theories, various modifications 

.£h have been found that produce, in particular, gauge/string duals with a mass gap, confinement, 

^ and supersymmetry breaking [3]. 

In this paper we continue our study of the condensates within gauge/string duality. In [5], we 
estimated the value of the gluon condensate with a result that is surprisingly close to the original 
phenomenological estimate of [I]. As known, QCD is a very rich theory supposed to describe the 
whole spectrum of strong interaction phenomena. The question naturally arises: How well does 
gauge/string duality describe other condensates? Here, we attempt to analytically determine 
a function Q and a parameter m\, which appear in nonlocal and mixed condensates, as an 
important step toward answering this question. 

The simplest quark-gluon mixed condensate is that associated with a dimension-5 operator 
constructed from the quark and gluon fields as qa^G^q, where a^ v is an antisymmetric combi- 
nation of 7-matrices. It is well known that it plays an important role in various QCD sum rules 
0. 



'Also at Landau Institute for Theoretical Physics, Moscow. 
1 For a review and list of references, see [2]. 
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The parameter m 2 , appears as a constant of proportionality in the conventional parametriza- 
tion [21 E] 

(gq^G^q) = m 2 (qq) , (1) 

where (qq) is a quark (chiral) condensate and g is a gauge coupling constant. 

One possible way to determine both the mixed condensate and the parameter TOq is to consider 
a nonperturbative gauge invariant correlator, also known as the nonlocal condensate^] 

ty(x 1 ,x 2 ) = {q(x 1 )U P (x 1 , x 2 )q(x 2 )) , (2) 

where Up(xi, x 2 ) is a path-ordered Wilson line defined as Up{x%, x 2 ) = Pexp[ig ds^- A^(x{s))\ 
Here s is a parameter of the path running from at x = x\ to 1 at x = x 2 . The path is taken 
to be a straight line. If one sets [6] 

*(xi,s 2 ) = {qq)Q(r), (3) 

then uiq is given by the coefficient of r 2 , with r = \x\ — x 2 \, in the expansion of the function Q 
as r — > 

Q(r) = l-^m 2 r 2 + O(r 4 ). (4) 

Note that this formula holds in Euclidean space. In Minkowski space it is modified by replacing 
r 2 -)• -r 2 [6]. 

Before proceeding to the detailed analysis, let us set the basic framework for the dual de- 
scription. We take the dual string spacetime as a product of five-dimensional space, with the 
Euclidean metric 

ds 2 = G nm dX n dX m = R 2 % (dt 2 + dx 2 + dz 2 ) , h = e BZ " (5) 

and some five-dimensional internal space X. Here s is a parameter whose value can be fixed 
from the heavy quark potential (Cornell model) or the slope of the Regge trajectory of p(n) 
mesons. We also take a constant dilaton and discard other background fields. In what follows, 
we assume a trivial dependence on the internal space X. Unlike other string duals, this model 
does share a few key features with QCD that singles it out and makes it very attractive for 
phenomenology: First, the model is a nearly conformal theory at UV, where it leads to the 
quadratic corrections [TJ [5]. Second, the model results in a phenomenologically satisfactory 
description of the confining potential [9l [10] . Finally, at finite chemical potential its extension 
provides the phenomenologically acceptable equation of state for cold quark matter jll| . 

2 The literature on the nonlocal condensates is vast. For more details, see, e.g., [B] and references therein. 
3 This is the ansatz. At this time we do not know equations that provide such a solution. So, we follow "the 
inverse scattering problem": first, we suggest a solution, then we look for its phenomenological relevance. 
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2 Calculating the Correlator 



To begin with, we will make an ansatz for computing the function Q within gauge/string duality 
whose justification, initially, is that it combines the ingredients at hand in the most natural way. 
Gradually, further evidence for the ansatz will emerge. 

Let us set the quark operators and the Wilson line on a four-manifold which is the boundary 
of a five-dimensional manifold. We will assume that the function Q is given in terms of the area 
(in string units) of a surface in the five-dimensional manifold by R 



Q = e~ s 

Here S is the area of a surface sketched in Figure 1. 



(6) 
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Figure 1: A surface in a five-dimensional manifold. The boundary is at z = 0. The surface is bounded 
by a curved profile of a static string stretched between the quark sources set at x = ±r/2 and a straight 
Wilson line along the x-axis. 



It is worth noting that a similar representation suggested in [12] has proven successful for 
studying the expectation value of the Polyakov loop in pure gauge theories. Recently, it was 
shown in [13] that it also reproduces the exponential decay of the correlator \£ at large separations, 
as expected in QCD. In the model with dynamical quarks it appeared in |14j . A difference is 
that here the Wilson loop goes along an internal direction on a five-dimensional boundary. 

2.1 Shape of Static String 

Following [9], we will now describe the shape of the static string stretched between the quark 
sources in the background geometry (|5|. First we set the quark operators at x = ±r/2 on the 
boundary (z = 0), as shown in Figure 1. Next, we make use of the Nambu-Goto action endowed 
with the background metric 



S= 2W / d2 ^V det Snmd a X n dpX m . 

In the static gauge £ = t and £ 2 = x, the action becomes 

(■r/2 



s = fr/ r dx^i + iz'f 

2VT J -(r/2) z l 



-(r/2) 

where q = R 2 /a'. A prime denotes a derivative with respect to x. 

4 We discuss some issues that arise in attempting to include a constant proportionality below. 
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Then, as in [9], we find the equation of motion for z 

zz" + 2(1- 5z 2 )(l + (z') 2 ) = 

as well as its first integral 

h 

const . 



(9) 



(10) 



The integration constant can be expressed in terms of the maximum value of z. On symmetry 
grounds, the function z{x) has a maximum at x = 0. It is z(0) = z , as shown in Figure 1. 

Since z{x) is an even function, we pick a fundamental domain defined by < x < r/2. 
Moreover, it is convenient to use the inverse function x+{z) subject to the boundary conditions 



1 



x+(0) = -r, x + (z ) = 0. 



(11) 



Using these boundary conditions, the solution to equation (10) in the fundamental domain is 

rj ,i 



x+(z) 



dnn 2 e A ( 1 -« 2 )(l-u 4 e 2A ( 1 -" 2 ))-( 1 / 2 ), 



(12) 



where A = sz 2 . Finally, the solution on the interval —r/2 < x < is simply x_(z) = — x + (z). 



If we set z = 0, then (12) becomes 

rj ,i 



duu\ x ^- u ^(i-uV^))-"/ 2 ), 



(13) 



which is the form in which it is written in [9]. 



2.2 Renormalized Area 



Having found the function x+(z) describing the string shape, we can now calculate the renor- 
malized area of the surface shown in Figure 1. 

To begin with, we choose the gauge £ — x and £ 2 = z. Then we substitute this into the 
action Q to obtain 



S 



e A« 2 

du — tt x + 
u z 



(14) 



V A Jo 

Here we have rescaled z such that z = z u. 

The integral is divergent at u = due to the factor z~ 2 in the metric ([5]). To proceed, 
therefore, we need to regularize it. A way to deal with this divergence is to cut off the integral 
at u = e. So, we have (up to terms vanishing for e — > 0) 



S 



TTZ 



r / 1 
2Ve 



1 + 



du 



e Xu2 [x, 



u- 



du 



r\ r 

2) + 2h * {V 



1 



(15) 



Now, by integration by parts, the action is 
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TT 



1 r 

2z e 



+ f 1 duu 2 ^ l - u2 \l - u 4 e 2X ^- u2 y [1/2 \V^XETfi(VXu) - u-h Xu2 ) 
Jo 



(16) 



where Erfi(u) is the imaginary error function. 

As in |13] , we use the modified minimal subtraction scheme to deal with the power divergence. 
So, we subtract ^-(^ — c), where c denotes a constant whose value must be specified from 



renormalization conditions. Finally, using ( 13 ) we get 
•l 



S 



IT 







duu 2 e x ^- u2 \l - n 4 e 2A ( 1 -" 2 ))- (1/2) (v / 7rAErfi(VAn) - u' 1 ^ + c) 



(17) 



2.3 Q and m 2 

Actually, the function Q = e~ s is written in parametric form given by equations (13) and (17). 
At this time it is not clear to us how to eliminate the parameter A and find Q as a function of 
r. We can, however, gain some important insights into the problem by considering two limiting 
cases. 

First, let us have a close look at r(A). According to [9j], it is a continuously growing function 
defined on the interval [0, 1]. The asymptotic behavior near zero is given by 



(18) 



where p = r 2 (|)/(27r) 3 / 2 . Since A — > means r — > 0, small A's correspond to small values of 



The asymptotic behavior near 1 is given by 



ln(l - A) + O(l) . 



(19) 



Thus, this region corresponds to large values of r. 

At this point a comment is in order. If we take A G [0, 1], then there exists the upper bound 



on the maximum value of z such that z < z c = l/\/s [9]. In view of the formula ( 19 ), this means 
that z — > z c as r — > oo. 

Once the behavior of r(A) is understood, we can use it to study the properties of the function 
Q at short and long distances. 

We begin with the case of small r. Expanding the right hand side of equation (17) in A, in 
next-to-leading order we have 



S 



IT 



(20) 



Then combining this with (18), we find the desired behavior of the function Q at short distances 

1 



Q = Qo[l-^ m o r +°( r 



where 



AcQ5p(p 2 -TT X ) 



(21) 

(22) 
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and 



Qc 



exp^il 



")}• 



(23) 



In a similar spirit, we can explore the long distance behavior of Q. It follows from (17) that 
in the neighborhood of A = 1 the renormalized area behaves as 



S 



_0_ 

"2vr 



(v^rErfi(l) - e + c) ln(l - A) + 0(1) . 



(24) 



Along with the relation (19), this means that the function Q decays exponentially at long dis- 
tances as 



Q 



-Mr 



where [3 



M 



^(^Erfi(l) 



2tt 



e + c 



(25) 



(26) 



This behavior is precisely analogous to what is expected in QCD. 

Having understood the two limiting cases, we are now able to make some estimates relevant 
to phenomenology. 

To make an estimate of the parameter m 2 ,, we need to fix a value for the constant c. First, 



we will fix it from the standard normalization of Q that is Q(0) = 1 [S]P] From (23), we find 
that c = vrp/2 ~ 1.3. Given the value of c, we can estimate the parameter m 2 ,, with the result 



ml 



1 



0.70 GeV 2 



(27) 



2gsp z 

Here we have used s ~ 0.45 GeV 2 and g ~ 0.94 as it follows from the fits to the slopes for the 
Regge trajectory of p{n) mesons [7] and the linear term of the Cornell potential [9j. 

There is a long history of attempts to estimate the value of m\ [2] . According to the original 
phenomenological estimate based on the QCD sum rules [15], which is widely accepted, it is given 
by m\ = 0.8 ±0.2 GeV 2 . Thus, our estimate is very satisfying at this point. Of course, there are 
other estimates which are also close to the original value. For comparison, lattice simulations 
|16| and the field correlator method [T7] both result in somewhat larger values, of order 1 GeV 2 Q 

To complete the picture, let us estimate the parameter M resulting from long distances. With 
our values for the parameters, we have 



M « 0.15 GeV. (28) 

This value is surprisingly closed to the pion mass that may be an indication that at long distances 
the correlator is dominated by the lightest meson contribution. 

On the other hand, our stringy construction suggests a natural normalization condition M = 
yfo, where a is the string tension (the coefficient of the linear term in the Cornell potential). 

5 The length scale £ = M^ 1 known as the correlation length has been also computed in |13] , 

6 While it may seem natural to assume this normalization in our framework, where it is equivalent to the fact 

that the surface shown in Figure 1 shrinks to a point as r — > 0, it requires a caveat. In section 3, we will propose 

a more refined way. 

7 At this point, we are somewhat formal and compare the results in different renormalization schemes. 
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For a typical value of a = 0.18 GeV 2 , it gives M ~ 0.42 GeV together with c ~ 4.0. This simple 
estimate allows us to draw the following conclusions: 

(i) The value of M is now close to the lowest energy level of the heavy-light mesons in the 
heavy quark effective theory [6j which is around 0.45 GeV. This suggests that at long distances 
the correlator is dominated by this lowest state. 

(ii) tuq is definitely larger than 0.8 GeV 2 (up to a factor of 2.6), and that it is of order 



2.1 GeV 2 . Interestingly, the lattice calculation of |18j and the instanton liquid model of [19] 
yield even larger values, about 2.4 — 2.5 GeV 2 . 

(iii) An additional free parameter is needed if one also wants Q(0) = 1. A simple way to 
introduce it is to use a slightly modified ansatz Q = e~ s /Qo, where Qq is the normalization 



constant defined by eq.(23). 



Some comments about the value of m 2 ,. The procedure for determining m 2 , is very sensitive 
to the procedure of extracting the power divergence. This is consistent with the literature, where 
the values show significant scatter (up to a factor of 3). One of the motivations for this work 
was to calculate it in a new nonperturbative approach that would shed light on this problem. 



Finally, let us present the results of numerical calculations. The parametric equations (13) 



and (17) predict a characteristic form for the function Q, as shown by the upper curve in Figure 
2. 




r(fm) 



Figure 2: Q as a function of r. The upper blue curve comes from the parametric equations (13) and 
\VJ\ . The lower red and green curves correspond to the models A and B, respectively. For all the cases, 
= 0.7 GeV 2 and M = 0.15 GeV. 



It is of great interest to compare this form with phenomenological models. Various possibili- 
ties of those have been discussed in the literature. Here we will consider two simple models, say A 
and B, for the vacuum distribution function / which in Euclidean space is related to Q by an in- 
tegral transform Q(r) = J °° dse~ sr / 4 /(s). Following [6], we take /a(s) = A exp{— M 2 / s — s 2 / a} 
and /s(s) = -Bexp{— M 2 /s — s/b}. Both these functions reproduce Q ~ e~ Mr for large r. The 
normalization is chosen so that the zeroth moment is ds f(s) = 1 or, equivalently, Q(0) = 1. 
In addition, the first moment must obey ds sf(s) = ttiq/A. We see that in the phenomenologi- 
cally important interval 0.1 fm < r < 1.5 fm our model is closer to the model B whose distribution 
function falls exponentially at large s. For r < 1 fm the agreement between the models is quite 
good. The maximum discrepancy occurred at r = lfm is of order 25 %. 
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3 Many Open Problems 



There is a large number of open problems associated with the circle of ideas explored in this 
paper. In this section we list a few. 

(i) Our analysis is not very accurate for various reasons. For one thing, the string breaking 
effect must be included in the calculation of the renormalized area S. What happens when a 
light qq pair is created? The surface is now modified by the string decay, like that of Figure 3. 
Thus, one might expect that this would require improvement of the analysis of section 2 at large 
distances. At this point it is worth noting that the string breaking scale is of order 1.1 — 1.2 fm. 
Therefore, it is tempting to say that the discrepancy with the phenomenological models of Figure 
2 for r > 1.2 fm is a result of the string breaking effect. 



Z 











-r/2 


r/2 



Figure 3: A surface modified by a light pair creation, there appears a cusp on the curve at x = 0. 



Also, the normalization of the nonperturbative correlator ^ (function Q) at r = is prob- 
lematic: one might worry that this is a territory of perturbative QCD. It is therefore interesting 
to see what happens if the correlator is normalized at a length scale ro 7^ 0. For doing so, let us 
set ro = 0.2 fm and see whether m\ has a reasonable value. As before, we take s ~ 0.45 GeV 2 



and g ~ 0.94. Using (13), we find numerically that a solution to r(A) = ro is A ~ 0.23. With 



this value of A, the renormalized area (17) vanishes at c ~ 1.13. Finally, we get from (22) 



ml « 0.60 GeV 2 , (29) 

with Q(0.2fm) = 1. This is still a satisfying and reasonable value. However, our analysis doesn't 
involve stringy quasiclassical corrections which may be important^] 
(ii) In the case of AdSs the parametric equations become trivial 

Note that S is independent of r, as required by conformal invariance. 

Now, a problem arises. If one tries to compute the function Q within the model based on a 
truncated AdS space, where the shape of a static string is the same as in AdS space until the 
string is long enough to reach the cutoff (IR brane location) [21], how can the expansion Q 
occur? □ 

8 Although we have no satisfactory explanation of why these corrections are small, it is interesting to note that 
the classical string computation of the heavy quark potential shows a remarkable agreement with the lattice data 

[3 [20]. 

9 In this model the computation of the gluon condensate along the lines of [5] is also puzzling. 
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(iii) An interesting observation of lattice simulations [22] is that at low temperatures (below 
T c ) the parameter ttiq is almost independent of the temperature. We can gain some understanding 
of this by writing the dual string spacetime either as ([5]), with t a periodic variable of period 
1/T, or as [23] 

ds 2 = R 2 ^(fdt 2 + dx 2 + r 1 dz 2 ) , h = e sz \ f = l-^y, (31) 

where z T = l/(wT). Since small r's correspond to small z 's, we expand the five-dimensional 
metric around the AdSs background by writing h and / as power series in z. The first term in 
the expansion Q is a constant which equals 1. It comes from the AdSs metric. The second term 
determines the parameter m,Q. It is due to the leading correction to the AdSs metric. In our case 
the correction is independent of the temperature, so the parameter m§ is independent, too. This 
reasoning, however, requires a caveat at high temperatures, where the expansion around AdSs 
is no more appropriate. 

(iv) In section 2 we cut off the integral over u. What if we had chosen to cut off the integral 
over z and then subtract |^ ( - — c) ? In that case the large r behavior of S remains unaffected 
since z —> 1/y/s as r — > oo, but z depends on r at smaller r so the small r behavior of S will be 
different [jf] Although we have no satisfactory explanation of why the regularization of section 2 
is more appropriate, we chose it because it leads to quite reasonable results. 

(v) As a final remark, we should point out that in holographic QCD mjj is treated as a 
free parameter whose value is fixed by phenomenology [23]. In such an approach one starts 
from a five-dimensional effective field theory action, somehow motivated by string theory but 
without higher derivatives terms (stringy a' corrections), and tries to fit it to QCD as much as 
possible. One might think of criticizing this approach on the grounds that it doesn't include (all) 
stringy a' corrections to the effective action (for example, see |25j). Therefore, it is tempting to 
determine the value of m,Q by saying that it can be done as a result of the a' corrections. As we 
have seen above, the effective string theory description provides evidence for this by showing the 
c/-dependence of Wq. 
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